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METHOD OF MONITORING BRAIN FUNCTION 
FIELD OF THE INVENTION 

5 The present invention relates to a method and system for the application of a mathematical 
model, and in particular a fixed order auto-regressive moving average model, to analyse 
electroencephalogram ("EEG") signals generated by a subject in order to assess and 
monitor the subject's brain function under conditions of health, disease and therapeutic 
io>frvention, 

10 

BACKGROUND OF THE INVENTION 

In clinical practice involving alterations in the level of consciousness, such as during the 
administration of sedatives or general anaesthetic agents, it is important to be able to 

15 quantify brain function. Most approaches rely upon the analysis of the brain's surface 
electrical activity, known as the electroencephalogram or EEG. In general the signal 
analysis method chosen is based on the statistical properties of the signal being analysed. 
The more closely matched the method used is to the signal properties, the more reliable, 
meaningful and accurate the resulting analysis will be. However these signal properties 

20 can only be known if the mechanisms and processes responsible for the generation of the 
signal are also known. 

To date none of these analysis methods of the brain's rhythmic electrical activity have 
incorporated any details of the underlying physiological mechanisms responsible for its 
25 genesis. Therefore their ability to measure, and thus monitor, brain function in the clinical 
setting is limited. 

This problem is overcome by the present invention which provides a more rational means 
of assessing and measuring brain function based on the detailed knowledge of the 
30 physiological mechanisms underlying the generation of the brain's surface rhythmic 
electrical activity. 
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The theory underlying the present invention considers the cortex of the brain as a single 
excitable spatial continuum of reciprocally connected excitatory arid inhibitory neurons 
interacting by way of short-ranged (intra-cortical) and long-range (cortico-cortical) 
connections. As such, the brain is seen as a dynamically evolving entity rather than a 
5 synthetic processing unit like a computer. 

Based on this theory, the characteristics of alpha rhythms arising as a consequence of the 
brain's neural connections can be closely represented by a mathematical model, and in 
particular, a fixed order auto-regressive moving average ("ARMA") model. The present 
10 invention derives specific values for the moving average ("MA") and auto-regressive 
4 ("AR") orders for the ARMA model based on the electrocortical transfer function. The 
electrocortical transfer function describes in a mathematical form the origin of the EEG 
readings taken of a subject. 

15 By applying EEG signals recorded from a subject to the fixed order ARMA model, 
coefficients can be obtained. To understand how these coefficients can be used to measure 
brain function, the equations defining the fixed order ARMA model are rewritten in the z- 
domain (complex domain) and are solved to obtain complex number solutions (called 
"poles") that are mapped onto the z-plane. These poles represent the state of the brain at 

20 the specific point in time when the EEG signal was recorded. Variations in the EEG 
signal, such as that induced by applying sedatives to the subject, can be detected as 
variations of the mean location of one or more poles on the z-plane. These variations can 
be interpreted to measure brain function or to indicate changes in the state of the brain. 

25 SUMMARY OF THE INVENTION 

According to the present invention there is provided a method for assessing brain state by 
analysing human electroencephalographic recordings using an eighth order autoregressive 
("AR") and fifth order moving average ("MA") discrete time model based on a theory of 
30 the underlying mechanism of generation of mammalian EEG activity. 
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The invention also provides a method for assessing brain state by analysing human 
electroencephalographic recordings using an eighth order autoregressive and fifth order 
moving average discrete time equation, taking a z-transform for said equations to obtain a 
z-domain equation, determining poles and zeroes in the solution of the z-domain equation 
5 and plotting the poles onto the complex plane. 



10 



15 



20 



The invention also provides a method of assessing the state of a mammalian brain 
including the steps of: 



(i) 

(ii) 

(iii) 

(iv) 



equation: 



obtaining an electroencephalogram ("EEG") from the brain; 
digitising the EEG to define a digitised data signal; 
segmenting the EEG into time frames of fixed length, y[n]; 
approximating each digitised time frame by a first equation: 



& 5 

y m = a k y^ n ~ *i + X b k U l n - *] 

*=0 



(v) solving the first equation to determine coefficients ai to ag and bo to b 5 ; 
, (vi) performing a z-transform on the first equation to obtain a second, z-domain 



-k 



Y(z) = U{z) = ^rrU(z) 



8 

i 



A(z) 



(vii) substituting each of the values of the coefficients into the z-domain 
equation; 

(viii) solving A(z) = 0 for z in the second equation to determine the poles; 

(ix) plotting the poles in the complex plane; 

25 (x) repeating steps (iv) to (ix) for each frame in the sample to determine 

clusters of poles in the complex plane; and 
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(xi) assessing the state of the brain by reference to the position and distribution 
of at least some of said clusters of poles as mapped in the complex plane. 

According to the present invention there is also provided a method of assessing the state of 
5 a mammalian brain including the steps of: 

(i) obtaining an electroencephalogram (EEG) from the brain; 

(ii) digitising the EEG to define a digitised data signal; 

(iii) segmenting the EEG into time frames of fixed length, y[n]; 

(iv) approximating each digitised time frame by a first equation: 



10 



O J 



(v) solving the first equation to determine coefficients ai to ag and bo to bs; 

(vi) performing a z-transform on the first equation to obtain a second, z-domain 
IS equation: 

Hz) = U(z) =-7^f (z) 

(vii) substituting each of the values of the coefficients into the z-domain 
equation; 

(viii) solving A(z) = 0 for z in the second equation to determine the poles; 
20 (ix) plotting the poles in the complex plane; 

(x) repeating steps (iv) to (ix) for each frame in the sample to determine 
clusters of poles in the complex plane; 

(xi) administering an intervention to the brain; 

(xii) repeating steps (i) to (x) at least once; 

25 (xiii) monitoring movement of at least some of said clusters of poles in the 

complex plane; and 
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(xiv). assessing the state of the brain by reference to movement of at least some of 
said clusters of poles as mapped in the complex plane. 

According to the present invention there is also provided a method of assessing the state of 
a mammalian brain including the steps of: 



(i) 

(ii) 

(iii) 

(iv) 



10 



obtaining a first electroencephalogram (EEG) from the brain; 
digitising the EEG to define a digitised data signal; 
segmenting the EEG into time frames of fixed length, y[n]; 
approximating each digitised time frame by a first equation: 

8 5 



Jfc=0 



15 



(v) 
(vi) 
equation: 



solving the first equation to determine coefficients ai to as and bo to b 5 ; 
performing a z-transform on the first equation to obtain .a second, z-domain 



5>» 



Y(Z) = 



_ k=0 



*=1 



A(z) 



20 



25 



(vii) substituting each of the values of the coefficients into the z-domain 
equation; 

(viii) solving A(z) = 0 for z in the second equation to determine the poles; 

(ix) plotting the poles in the complex plane; 

(x) repeating steps (iv) to (ix) for each frame in the sample to determine 
clusters of poles in the complex plane; 

(xi) obtaining a second EEG from said brain at a later time; 

(xii) repeating steps (ii) to (x) in relation to the second EEG at least once; 

(xiii) monitoring the movement of at least some corresponding clusters of poles 
in the complex plane derived from the first and second EEGs respectively; and 
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(xiv) assessing the state of the brain by reference to movement of at least some of 
said clusters of poles as mapped in the complex plane. 

For the methods above, an EEG may be obtained and recorded before it is processed. The 
5 recorded EEG can therefore be processed at any time after it has been recorded or it can be 
used as a reference for comparisons with other EEGs at a future point in time. 
Alternatively, an EEG may be obtained and processed on-the-fly such that an EEG is 
repeatedly obtained over consecutive and constant time intervals, and where each time 
interval may overlap with the immediately preceding time interval. The EEG obtained for 
10 each time interval is immediately processed by the methods described above. 

Preferably, step (x) is repeated up to 100 times so that there are a plurality of poles in each 
of the said clusters. Also step (x) may be repeated continuously to track the motion of the 
poles from each segment. 

15 

Preferably further, the method includes a step of taking the centroid of the poles for each 
cluster of poles, and monitoring and comparing the movement of the centroids. 

The present invention further provides a system for performing the above methods. The 
20 present invention further provides computer readable media having computer program 
instructions stored thereon which, when executed by a computer, perform the methods 
described above. 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 

Preferred embodiments of the present invention are hereinafter described, by way of 
example only, with reference to the accompanying drawings, wherein: 

FIG. 1 is a table that defines the range of all the theoretical parameters (i.e. the 
numerical values of all anatomical and physiological parameters) that are used by some of 
30 the equations referred to in this specification to generate parameter sets that give rise to 
stable physiological alpha activity; 
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HG. 2 is a diagram showing an example of the 8 poles derived from one segment 
of an EEG signal plotted onto the z-plane; 

FIG. 3 is a diagram showing an example of the cumulative positions of the 8 poles 
derived from several segments of an EEG signal plotted onto the same z-plane. 
5 FIG. 4 is a diagram showing the schematic representation on the z-plane of the 

predicted effects of increasing the strength of neuronal population inhibitory-inhibitory 
and inhibitory-^excitatory synaptic interactions; 

FIG. 5 is a diagram showing the schematic representation on the s-plane (which can 
also be referred to as the Laplace or Fourier plane) of the predicted effects of increasing 
10 the strength of neuronal population inhibitory-inhibitory and inhibitory-+excitatory 
synaptic interactions; 

FIG. 6 is an example of the upper right quadrant of the z-plane in a pole-zero plot 
for a typical subject before (as shown in black) and after (as shown in grey) the 
administration of the benzodiazepine, alprazolam, on the subject; 
15 FIG. 7 is an example of a detailed view of a region of the upper right quadrant of 

the z-plane between 8 to 13 Hz of a pole-zero plot for a typical subject before (as shown in 
black) and after (as shown in grey) the administration of the benzodiazepine, alprazolam, 
on the subject. 

FIG. 8 is an example of the upper right quadrant of the z-plane in a pole-zero plot 
20 for a typical subject before (as shown in black) and after (as shown in grey) the 
administration of a placebo on the subject; 

FIG. 9 is an example of a detailed view of a region of the upper right quadrant of 
the z-plane between 8 to 13 Hz of a pole-zero plot for a typical subject before (as shown in 
black) and after (as shown in grey) the administration of a placebo on the subject. 
25 FIG. 10 is a diagram showing one example of the physical setup of components in 

the system; 

FIG. 11 is an example of an EEG signal typically recorded from a subject over the 
course of a 10-second interval. 
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DET AILED DESCRIPTION OF PREFERRED EMBODIMENTS OF THE 
INVENTION 

The alpha rhythm is arguably the most obvious recordable feature of the intact human 
5 brain. While the exact basis for its genesis is still controversial it is widely believed that it 
arises as a consequence of one or more of the following mechanisms: 

• endogenous or exogenous (thalamic) pacing of cortical neurons; 

• oscillatory activity generated through the reciprocal interactions of excitatory 
(pyramidal) and inhibitory (interneuron) cortical neuronal populations; or 

10 • boundary dependent standing wave generation (Schumann like resonance) due to 
long range cortico-cortical connectivity. 

However none of these mechanisms are sufficient, either separately or taken together, in 
explaining the physiological genesis of the alpha rhythm. 

15 

A theory of alpha electrorhymogenesis (as discussed in Liley et al. Network: Comput. 
Neural Syst. 13 (2002) 67-113, the contents of which are hereby incorporated in this 
specification) is based upon a detailed spatially continuous two-dimensional mean field 
theory of electrocortical activity. According to this theory, the brain acts as a white noise 
20 filter to its electrical neural input and the alpha rhythm arises as a result of the filtering of 
input signals going to the cortex. The filter properties are determined by the bulk 
(macroscopic/large-scale) anatomical and physiological properties of excitatory and 
inhibitory cortical neurons. 

25 In this theory, inhibition is conceived as having an important role in determining the 
properties of the "cortical filter" and thus the spectra of the alpha rhythm generated. In 
particular the selective modification of the strength of cortical inhibitory action by 
benzodiazepines, such as alprazolam, is associated with specific changes in the properties 
of this filter. As such, it is found that the strength and form of the population 

30 inhibitory-inhibitory synaptic interactions are the most sensitive determinates of the 
frequency and damping of the emergent alpha band oscillatory activity. Such behaviour 
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10 



arises principally because local inhibitory-inhibitory and local inhibitory-excitatory loop 
delays that are associated with physiologically and eiectroencephalographically plausible 
alpha activity are longer than the corresponding local (intra-cortical) and long-range 
(cortico-cortical) excitatory-^excitatory loop delays. 

This theory differs from other macroscopic continuum theories in that the time course of 
the unitary inhibitory post-synaptic potential ("IPSP") is described by a third order 
differential equation. Lower orders are theoretically found to be unable to support any 
appreciable or widespread alpha band activity. 



The principal state variables modelled under this theory are the mean soma membrane 
potentials of local cortical populations of excitatory and inhibitory neurons. The local field 
potential, and hence the EEG or electrocorticogram ("ECoG") signal, is regarded as being 
linearly related to the mean soma membrane potential of the excitatory neurons. This 

15 theory can be cast as a set of coupled non-linear one-dimensional partial differential 
equations that incorporate the major bulk anatomical and physiological features of cortical 
neurons and includes cable delays, neurotransmitter kinetics and cortico-cortical and intra- 
cortical connectivities. The spontaneous alpha rhythm is theorized to predominantly arise 
as a consequence of the local linear properties of the cortex. For this reason in the current 

20 formulation spatial effects have been restricted to one dimension. 

In accordance with this theory the following non-linear equations (Equations 1, 2, 3 and 4) 
mathematically represent the brain's electrical activity, as further described in Liley et al. 
Network: Comput. Neural Syst. 13 (2002) 67-113: 

r aM* 1 t) = _ *(»»*) + •-<*>i P .(«»*) + Equation 1 
at 

{§t +% ) Ie ^ X,t) = r «?«^( 1 H^.(ft.) + #x,t)+p.(x,t)> Equation2 

(j^ + 7< ) I { (x,t) = r^i&qpilHNfSM+ptix^)} Equation 3 

^ (i^ + vaJ 4{x,t) - v 2 ^^^- = vAN" (vA + I^S e {h e ) Equation 4 
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where 



N f = ( N f e ,Ng) T , N°=(NZ,Nfl, #«fo.*) r . A = ^ag(A w ,A d ), r = dtog(r e ,r l ), 
^W = ^(^(^^(M). P.=(p«.pJ. p t =(p„P„) r and /is the identity matrix, 



with: 



5 



5,(hj) = 57"(1 + exp[-V^(^ - 1 Equation 5 

*i(M = (fcr-M/W-^l Equation 6 

where/ / = e, i. Figure 1 is a table that defines the range of all the theoretical parameters 
(i.e. the numerical values of all anatomical and physiological parameters) that are used by 
10 the above equations to generate parameter sets that give rise to stable physiological alpha 
activity. The ranges in Figure 1 refer to the intervals from which uniform parameter 
deviates were generated. 

These non-linear equations need to be transformed into their linear equivalent in order to 
15 be solved. To determine theoretically whether the alpha rhythm can be understood in 
terms of a white noise fluctuation spectrum the above equations are linearized about 
spatially homogeneous singular points. For a given set of parameters these singular points 
can be obtained by setting all spatial and temporal derivatives to zero and solving for h e . 
In general these singular points, h' e , are solutions to the following equation: 



20 



25 



F(h e (q), q) = 0 Equation 7 

where q represents a vector of model parameters and F(.) is obtained from Equations 1, 2, 
3 and 4. 

Linearizing Equations 1, 2, 3 and 4 about the spatially homogenous singular point h\ and 
transforming to the Fourier domain yields the following equation: 
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r e D{k 7 u:q) 
= G e (k 9 v:q)P(k,uj) Equation 9 

where k and (o are wave number and angular frequency respectively. Loosely speaking, k 
specifies the reciprocal of the characteristic physical scale over which oscillations of 
frequency <o occur. H e (k,<o) is the Fourier transform of the mean soma membrane potential 
of excitatory neurons h e (x,t). he(x t t) has been shown to be proportional to the surface 
recorded electrical activity, the EEG, of the brain. The function G e is the electrocortical 
transfer function, q is a vector of parameters and P(k,(o) represents the spatio-temporal 
form of cortical input. 

The terms N(k,G>,q) and D(k,a),q) in Equation 8 can be expressed as the following 
Equations 10 and 11, where the corresponding parameters for the parameter vector V in 
Equation 8 has now been explicitly identified in Equations 10 and 11. 

N(k,u>) = {(iu; + 7<) 2 (*" + W^-^ Equation 10 

-w cc Qc(N: c k(k + n>M + + + k 2 })] 

-w^tiMcNPQeQiiN^KiK + iu/v) + N^{(k + iu/v) 2 + k 2 }) Equation 1 1 

w fJ = exp(l)r^^^(h;)/(r,v77 i ) Equation 12 
Tfe = l+exp(l)r e (^ 

^Ml^iNfjSdhm^hr - h'f«\) Equation 13 

Qi = ^S s /dhj\ hjssah j Equation 14 

From Equation 10, the highest order in co is 5, which corresponds to the moving average 
order of the ARMA model. From Equation 11, the highest order in m is 8, which 
corresponds to the auto-regressive order of the ARMA model. 



Equations 8, 10 and 11 can be rewritten in a summary form as Equations 12, 13 and 14, as 
shown above. 
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Equation 14 can be rewritten as a difference equation, as shown in Equation 15, which 
represents a linear time invariant discrete time system: 

8 5 
yW = ~£ yi n ~ *] + 2j b k U ^ n ~ k 1 Equation 15 

k=l k=0 

5 

where y[n] is the digitised EEG signal, u[n] is a Gaussian white noise process and a k and 
b k are coefficients to be determined for a given EEG time series. 

Equation 15 represents an (8,5) order ARMA model, where the specific values of the 
10 orders are derived from Equations 10 and 1 1. The 14 coefficients from the ARMA model 
can be determined using any of the large number of commercially or freely available 
ARMA software modelling packages, such as the ARMASA Matlab Toolbox software by 
P.M.T Broersen (Delft University of Technology). 

15 To understand how the 14 coefficients so obtained can be used to measure brain function, 
Equation 15 is rewritten in the z-domain by taking the z-transform. Thus Equation 15 can 
be equivalently written in the z-domain as: 



20 



Y( Z ) = -^g U(Z) = ^riU(z) Equation 16 



Solutions to A(z) - 0 in Equation 16 will give the system poles and solutions to B(z) = 0 in 
Equation 16 will give the system zeros. In general, these solutions are complex with 
\z\ < 1 . The maximum power of the denominator in Equation 16 suggests that there are 8 
unique system poles. The eight complex solutions to A(z) = 0 are then plotted on the z- 
25 plane. 
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The location of the eight poles derived ftpm Equation 16 represent the state of the brain as 
determined by the EEG signals recorded over the particular time interval. With reference 
to Figure 2, the set of eight poles from one EEG segment 20 is plotted on the z-plane 21. 
With reference to Figure 3, the application of Equation 15 to subsequent segments of EEG 
5 from the entire EEG sample allows further sets of eight poles 31 to be determined and 
plotted onto the same z-plane 30 that contains the poles plotted from the preceding EEG 
samples 31. 

. Theoretically, the results in Figure 3 can be interpreted with reference to Figure 4. Figure 
10 4 is the schematic representation of the predicted effects of increasing the strength of 
neuronal population inhibitory-inhibitory and inhibitory-*excitatory synaptic interactions. 
The filled circle 40 approximately represents the theoretical loci of the dominant poles 
associated with electroencephalographically plausible eyes-closed alpha activity. In other 
words, the filled circle 40 represents the most weakly damped pole, which can be thought 
15 of as corresponding to the most dominant oscillatory component making up the human 
alpha rhythm. For any recording of human alpha rhythm, a frequency analysis using a 
Fourier transform would reveal the approximate frequency of this dominant pole or the 
dominant oscillatory component. The arrows 41 and 42 in Figure 4 indicate the mean 
predicted direction of the motion of these poles in response to increases in 

20 inhibitory-inhibitory (3a/ /dN$ ) 41 and inhibitory-*excitatory (da* IdN* ) 42 synaptic 

strength. 

Figure 5 shows the same information as Figure 4 as plotted on the s-plane (Laplace or 
Fourier plane) rather than the z-plane. The filled circle 50 approximately represents the 
25 theoretical loci of the dominant poles associated with electroencephalographically 
plausible eyes-closed alpha activity. The arrows 51 and 52 in Figure 5 indicate the mean 
predicted direction of the motion of these poles in response to increases in 

inhibitory-inhibitory (da)* /&V,f ) 51 and inhibitory->excitatory {dm* ldN? e ) 52 synaptic 
strength. 

30 
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The decay rate is related to the sharpness of the resonance of the dominant oscillatory 
component in the recorded EEG signal. Increasing decay rates would correspond to the 
broadening of the alpha resonance in human EEG recordings. In Figure 5, the arrow 51 
represents an increasing decay rate as it moves away from the filled circle 50, which 
5 represents a single pole. The arrow 52 represents a decreasing decay rate as it moves away 
from the pole 50. In Figure 4, the closer that the pole 40 moves to the boundary of the unit 
circle 43 the smaller the decay, whereas the further the pole moves from the boundary 43 
the larger the decay. Referring to Figure 4, an anti-clockwise motion of the pole 40" 
implies that the pole's frequency increases, whereas a clockwise motion of the pole implies 
10 that the frequency decreases. Thus, according to this theory of alpha rhythm generation, 
changes in the decay rate(s) associated with the dominant oscillatory component gives 
information in addition to that which can be obtained using the Fourier analysis of human 



15 Figures 6 and 7 show the position of the poles and zeros in a pole-zero plot for a typical 
subject before (as shown by "-BZ" in black) and after (as shown by "+BZ" in grey) the 
administration of the benzodiazepine as an oral dose of alprazolam. A similar pole-zero 
plot for a typical subject is shown in Figures 8 and 9 before (as shown by "-PL" in black) 
and after (as shown by "+PL" in grey) the administration of a placebo. Where they apply 

20 in Figures 6, 7, 8 and 9, poles are indicated with and zeros indicated with "O". 

Figures 7 and 9 show in more detail the region of the z-plane, from Figures 6 and 8 
respectively, corresponding to a region of 8-13 Hz activity. The important features to note 
are: 

25 (i) the distinct groupings of poles and zeros populating the z-plane; 



EEG. 



30 



(ii) 
(iii) 



distinct populations of poles having frequencies lying in the range 8-13 Hz; 
clear differences between the location of the centroids of the alpha poles 
before .(indicated as "-BZ" in Figures 6 and 7, with and "O" symbols 
marked in black) and following (indicated as "+BZ" in Figures 6 and 7, 
with and "O" symbols marked in grey) administration of the 
benzodiazepine; and 
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(iv) insubstantial differences between the location of the alpha poles before 
(indicated as "-PL" in Figures 8 and .9, with "+"and "O" symbols marked 
in black) and following (indicated as "+PL" in Figures 8 and 9, with 
and "O" symbols marked in grey) administration of the placebo. 

5 

It is evident from Figure 7 that the variability of the alpha pole location after the ingestion 
of the benzodiazepine is more pronounced that in the other three conditions. Compared to 
the placebo condition, as shown in Figure 9, it can be seen that alprazolam causes a 
significant shift in the most weakly damped pole constituting the alpha rhythm, such that 
. 10 its corresponding frequency and damping both increased. In the absence of any other poles 
this implies that alprazolam causes the alpha spectrum to shift to the right and broaden, as 
shown in Figures 7 and 9. Variations in the mean location of these clusters of poles and 
their distribution will result as a consequence of ongoing normal behaviour, disease 
progression and state, or pharmacological or therapeutic interventions or manipulations. 
15 As such, the mean location of the poles is a measure of brain state and the movement of the 
mean location of these poles is a measure of changes in brain state or function. 



Figure 10 is a diagram showing a preferred embodiment of the physical components of the 
system. The EEG signals are detected by multiple electrodes 101 from the scalp of a 
20 subject 102, where the electrodes are positioned preferably arranged according to the 
international 10:20 standard system with the addition of mid-point electrodes as necessary. 
Preferably, the EEG signals are recorded referenced to linked ears using 64 scalp 
electrodes attached to an electrode cap using the nasion as a ground. Other techniques for 
recording the EEG can also be used. 

25 

It is preferable that the EEG is analysed as a sequence of overlapping fixed length 
segments. This technique is further described with reference to Figure 11. Figure 11 
shows an example of a typical EEG signal 201 recorded from a subject over a period of 10 
seconds. To illustrate the preferred sampling technique, it is assumed that recorded EEG 
30 signals are digitised and segmented into fixed, overlapping, 2 second segments. The first 
EEG segment 202 extends from 0 to 2 seconds. The next EEG segment 203 overlaps with 
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the preceding EEG segment (e.g. from 1 to 3 seconds). This process is repeated for all 
subsequent segments of the EEG signal. This example is based on the assumption that the 
degree of signal overlap for all EEG segments is 50% of the immediately preceding 
sample. While this only reflects the best segmenting practice, it is possible to segment the 
5 EEG signals into any constant time interval and with any degree of overlap. The EEG 
signal segment is typically sampled (i.e. digitised) at anywhere between 200 to 500 
samples per second. 

Referring to Figure 10, the EEG signal is transmitted as an analog signal via multiple 

10 independent electrical connections 103 that connect each electrode 101 to the input ports 
104 of an analog filter and amplification device 105. The signals from each electrode, 
when grouped together, constitute the aggregate EEG signal which is then amplified and 
filtered in the analog filter and amplification device 105. The analog EEG signal is sent to 
an analog/digital ("A/D") converter 106, which digitises the filtered analog EEG signal. 

15 Preferably, filtering of the digital EEG signal is performed by the analog filter and 
amplification device 105, which removes any 50 Hz artefacts and other sources of noise 
that may contaminate subsequent signal analysis. However digital filtering can also be 
performed by software running on the central processing unit ("CPU") of a personal 
computer ("PC"), The digitised EEG signal is sent to a PC 107 vtea data connection 108, 

20 which includes a serial or parallel connection. Upon entering the PC, the digitised EEG 
signal is se ? nt to the CPU 109 via internal data bus connections. The CPU 109 controls a 
memory module 110, which may comprise of random access memory ("RAM") 
components for short-term or temporary storage and recall of data, or a hard disk or 
another device that provides more permanent storage. The software for processing the 

25 digital EEG signal is also stored in either the RAM or hard disk of the memory module. 

The system may process the digitised EEG signal on-the-fly, such that an EEG is repeated 
obtained over consecutive and constant time intervals, and where each time interval may 
overlap with the immediately .preceding time interval as described above. The EEG 
30 obtained for each time interval may be temporarily stored in the RAM memory 
components of the system before it is processed shortly after it has been put in the RAM 
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and removed from the RAM after the EEG for that time interval has been processed. 
• There may be more than one EEG stored in the RAM at any time, which corresponds to 
the EEGs obtained for different time intervals. 

5 The digitised EEG signal may also be obtained and recorded before it is processed. The 
digitised EEG signal may be recorded on more permanent forms of storage, such as a hard 
disk, tape drive or a compact disc ("CD"). The recorded EEG can therefore be processed 
at any time after it has been recorded or can be used as a reference for comparisons with 
otitpr EEGs (that may be recorded from the same subject or also from different subjects) at 
10 a future point in time. 

Referring to Figure 10, the CPU 109 runs software to perform ARMA modelling on the 
digital EEG signal and to calculate the 14 ARMA coefficients according to Equation 15. 
This may be done using the "ARMASA Matlab Toolbox" software by P.M.T Broersen 
15 (Delft University of Technology), or any of the large number of commercially or freely 
available ARMA software modelling packages. 

Upon determining these 14 coefficients, the CPU 109 uses software, which may be the 
same software package as described above, to calculate and graphically plot the 8 pole 

20 positions on the z-plane for that EEG segment. The software instructs the CPU 109 to 
send the graphical data generated by the software to a display device 111 controlled by the 
CPU 109, in which the display device 111 may be connected to the CPU 109 via internal 
data bus connections. The display device 111 generates a visual representation of the 
information within the graphical data generated by the software, which may be in the form 

25 of a graph or chart as shown in Figures 6, 7 8 and 9. 

Although Figure 10 shows that the system may be implemented with the assistance of 
additional hardware components, it is possible to implement some of the features provided 
by the hardware using software. For example, with reference to Figure 10, the functions 
30 provided by the analog filter and amplifier 105 and A/D converter 106 can be implemented 
using software. As such, it is possible to implement the process of filtering and amplifying 
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an analog EEG signal, A/D conversion (digitising the EEG signal), segmentation of the 
digitised EEG signal, storage of the digitised EEG signal, solving of the 14 ARMA 
coefficients for each segment of the digitised EEG signal and generating a graphical plot or 
other visual representation of a single set or consecutive sets of 8 poles derived from one 
5 or multiple segment of the digitised EEG signal respectively, as one software package. 

The reference to any prior art in this specification is not and should not be taken as an 
acknowledgement or any form of suggestion that prior art forms part of the common 
general knowledge in Australia. 

10 

Many modifications will be apparent to those skilled in the art without departing from the 
spirit and scope of the invention. 

15 Dated this 20 th day of January, 2003 

Swinburne University of Technology 

by DA VIES COLLISON CAVE 
20 Patent Attorneys for the Applicant 
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